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Abstract 

The boundary conditions corresponding to the Matsubara formalism 
for the T > partition function may be introduced as constraints in 
the path integral for the vacuum amplitude. We implement those con- 
straints with time-independent Lagrange multipliers and, by integrat- 
ing out the original fields, we obtain an alternative representation for 
the partition function, in terms of the Lagrange multipliers as dynami- 
cal fields. The resulting functional integral has the appealing property 
of involving only (i-dimensional, time independent fields, and looks like 
a nonlocal version of the classical partition function. We develop this 
formalism within the context of the scalar and Dirac fields. 

Many years ago [T], the theoretical foundations for a systematic treatment 
of Quantum Field Theory (QFT) at finite temperature (T > 0) and its 
relativistic generalization [2] were laid down. That approach, the now called 
Matsubara (or imaginary time) formalism has been very successful indeed 
in allowing for the evaluation of thermal effects in QFT, both in the High 
Energy [3] and Condensed Matter Physics realms. Besides, the notion of 
Abelian and non Abelian gauge fields, with all its consequences for particle 
physics [31 SI [5] could also be set up in the T > context in a quite natural 
way. 
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A fundamental property introduced by this formalism is the imaginary- 
time periodicity (antiperiodicity) conditions for the bosonic (fermionic) fields, 
something that may be already discovered at the level of the partition func- 
tion, Z(j3), for a system at a temperature T = 1/(3: 



Assuming, for the sake of simplicity, that the system has only one degree 
of freedom, q(r), the expression above may be written more explicitly as 
follows: 



where we have introduced the rotating basis \q,r), for a fictitious imaginary 
time evolution. Then the standard path integral construction for the transi- 
tion amplitude between the times r = and r = (5 may be applied, to write 
the partition function in the Matsubara formalism: 



where one only considers phase-space paths q(r), p(r) (r e [0,/?]) such that 
q(0) = q{(3). There is no need to impose any boundary condition on the p(r) 
paths 0. 

In the usual case, one may integrate out the canonical momentum vari- 
ables, what leads, for the case of a d+ 1-dimensional field theory, to a model 
defined on S 1 x R d 7 where the radius of S 1 is proportional to the inverse 
temperature (3. In Fourier space, the corresponding frequencies are the usual 
discrete Matsubara frequencies. 

A characteristic feature of the Matsubara formalism (shared by the real- 
time formulation) is that the introduction of a time dependence for the fields 
is seemingly unavoidable, even when one limits oneself to the calculation of 
static thermal properties. 

With the aim of constructing a new representation where only static fields 
are involved, in this letter, we introduce an alternative way of dealing with 
the T > QFT partition function. It is inspired by a recent paper in which 
a constrained functional integral approach is used to implement the effect 
of fluctuating boundaries in the Casimir effect [6]. In the present context, 
this allows one to introduce the periodicity conditions by means of Lagrange 
multipliers (c/-dimensional when the field lives in d + 1 dimensions). Then 

1 Albeit a more symmetric form for the boundary conditions in the path integral for 
Z (f3) can be used [7]. 



Z(P) = Tr[e~ pH ] . 



(1) 






(3) 
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the original fields can be integrated, what leaves a functional depending only 
on the ci-dimensional Lagrange multipliers. 

We shall introduce the new method by considering first the case of a single 
degree of freedom, extending it to more complex systems afterwards. We 
begin by writing the phase-space path integral for Z , the vacuum functional 
for a single quantum degree of freedom: 



Z = J VpVq exp{- S [q(r) , p(r)]} 



(4) 



where S is the first-order action, So = f_°° dr £, with £ = —ipq + H(p, q), 
with H denoting the Hamiltonian, assumed (for the sake of simplicity) to be 
of the form: H(p, q) = T{p) + V(q). 

Of course, Z may be regarded as the limit of an imaginary-time transition 
amplitude, 

Z = lim (q Q ,T\q ,-T) 

= lim y \(q \n)\ 2 e- 2TE - = lim | (q Q \Q) \ 2 e~ 2TE » (5) 

n 

where we have introduced \n), the eigenstates of H, H\n) = E n \n), and q , 
the asymptotic value for q at T — >• ±00 (usually, q = 0). E is the energy 
of |0), the ground state. 

We have already recalled the usual functional integral for the thermal 
partition function Zq{(5) in the Matsubara formalism. Let us now show how 
one can write an alternative expression for the same object, starting from the 
vacuum transition amplitude, Z , and imposing a periodicity constraint on 
the paths. We first use the superposition principle, introducing decomposi- 
tions of the identity at r = and r = f3, in order to write Z in an equivalent 
way: 

Z = lim / dq 2 dqi{q ,T\q 2 ,f3) (q 2 ,f3\qi,0) (qi,0\q o ,-T) , (6) 

T-*oo J 

or, in a path integral representation, 

rq(T)=qo 



Zq = lim / dq 2 dq\ / Dp V>q e J P T 

T ^°° J Jq(P)=q2 

rq{P)=q2 rq(o)=qi 

x / VpVqe-^ drC / VpVqe-$-T dTC . (7) 

Jq(0)= qi J q (-T)=q 

The representation above is useful in order to show that the naive procedure 
of introducing a periodicity constraint for just the coordinate q(r) is not the 
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right procedure to obtain Z(f3) from the path integral for Z . Indeed, we 
can see that 

J VpVq 5{q{(3) - <z(0)) exp { - S [q(r), p(r)} } 

lim e - £ o(2T-/3) f dqimqi) m qi ^\ qh o) , ( 8 ) 



and taking the ratio with the (unconstrained) vacuum functional, 

= e p / dgiK0|gi)| (gi,/%i,0) , (9) 



f VpVqe' 



So 



where it is evident that we cannot extract a Z{(3) factor. There is, however, 
a simple alternative that does the job: to reproduce Z (j3) we impose period- 
icity conditions for both phase space variables, i.e., we introduce the object, 
Z s (0) defined as follows: 

z (g) I S(q(P) - q(0)) 5(p((3) - p(Q)) e~ s ° 

s[ ^ ) J VpVqe~ s ° ' ^ ' 

Using the superposition principle, we can see that: 



VpVq 5{q(p) - q(0))6(p{P) - p(0)) e" 5 ° = lim / dp x dq x 
(qo,T\ Pl ,P) (pi,P\qi,/3) (gi,/%i,0) (q h 0\ Pl ,0) ( Pl , 0\q , -T)l (11) 



or 

I VpVq5(q(p)-q(0))5(p(p)-p(0))e~ So = lim e - ^ 2 ^ 



x I ^ 1 (qo\0)(0\p 1 )(q 1 ,P\q 1 ,0)(p 1 \0)(0\q ) 



= hm e~ E ^ T -V |(g |0)| 2 f d Ql (ft,/3| ft ,0) . (12) 
We have used the property 

(pi, (gi,0|pi,0> = (pi|gi)(gibi) = ^ . (13) 
Then we conclude that 

Z s (f3) = Tr[e-e-- 6: ] (14) 
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where : H : denotes the normal-ordered Hamiltonian operator, i.e. 



:H:= H — Eq. (15) 

The conclusion of the previous derivation is then that, by imposing peri- 
odicity in both phase space variables, and discarding /^-independent factors 
(as required by the normalization) we obtain Z s (/3), the partition function 
corresponding to the original Hamiltonian with energies redefined to for 
the ground state. The subtraction of the vacuum energy is usually irrele- 
vant (except for some exceptional situations), as it is wiped out when taking 
derivatives of the free energy to calculate physical quantities. 

As a test, let us apply the previous results to the case of the harmonic 
oscillator in one spatial dimension. To proceed, we exponentiate the two 
5-functions, by means of Lagrange multipliers £1 and £2- They are just real 
variables in this case, and we are thus lead to a new expression for Z s (@), 
which may be written as follows: 



2vr 



?W6 (s(t-/3)-S(t)) +p(r)6 (<5(t-/3)-5(t)) 



i f +°° dr 

J OO 

x e 



where the constant M is used to denote the (Zq)" 1 factor. 

The integral over p and q (which we shall perform first) is obviously 
Gaussian. Indeed, interchanging the order of integration, we see that (jj]) 
may be rewritten as follows: 

X 1^,°° dr Qa(T)K. ab Qa{T)+i drj a {r)Q a {r) 

where we have introduced Q = (Q a ) (a = 1, 2) for the canonical coordinates, 
such that Qi = q and Q 2 = p. Also, j a ( T ) = (5(r — (3) — S(t)), and )C a b 
are the elements of the 2x2 operator matrix AT, given by: 

The result of integrating out Q a , may be written as follows: 



Z s {(3) = A/-27r(det/Cp J ^L e ~^ M ^ 



(19) 
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with 

M = 2ft(0) - - Q(-p) , (20) 

where fi(r) denotes the inverse of the operator /C of ([TBI) : its explicit form 
may be easily found to be: 

fsgn(r) \ W , T | 



(sgn = sign function). Equation ([2~Tj) can be used in (l20l . to write: 

M = ( W ~ 1 ° ) (n B (o;) + I)" 1 , (22) 

where 

n B (cu) = (e^ - I)" 1 (23) 

is the Bose-Einstein distribution function (with the zero of energy set at the 
ground state). 

The resulting integral over the £ a variables becomes: 

Z,(p) = —e 2 I"jM+'] . (24) 

J 27T 

This integral is over two real variables £ , which are 0-dimensional fields, 
one dimension less than the + 1 dimensional original theory. This inte- 
gral may be compared with the one corresponding to classical statistical 
mechanics. To that end, we evaluate the partition function in the classical 
(high-temperature) limit. In that limit, we approximate the integrand and 
Z s {j3) becomes: 



Zs{(5) ~ J S e "^ (?1 ' 6) (/9<<1) 



(25) 



where: 

= + (26) 



We see that (1251) corresponds exactly to the classical partition function for a 
harmonic oscillator, when the identifications: £i = p (classical momentum), 
and £2 = <? (classical coordinate) are made 

ZM ~ [^e-efa+W) (P«l). (27) 

./ 27T 
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On the other hand, had the exact form of the integral been kept (no ap- 
proximation), we could still have written an expression similar to the classical 
partition function, albeit with an 'effective Hamiltonian' H e ff(£i, £2): 




/3#e//(«l,£2) 



(28) 



where: 



(29) 



This shows that the quantum partition function may also be written as a 
classical one, by using a /3-dependent Hamiltonian, which of course tends to 
its classical counterpart in the high-temperature limit. This representation 
is also valid for interacting theories, see below. 

By integrating out the Lagrange multipliers in the (exact) expression for 
the partition function ( |24l) . we obtain: 



which is the correct result. In what follows, to simplify the notation, we 
shall omit writing the V subscript in Z{(3), assuming implicitly that one is 
dealing with the normal-ordered Hamiltonian. 

An important fact that has emerged from an analysis of the classical 
(high-temperature) limit: the Lagrange multipliers have a physical interpre- 
tation. The multiplier associated to the periodicity condition for q plays the 
role of a classical momentum, while the one corresponding to the periodic- 
ity for the momentum becomes a generalization of the classical coordinate. 
The same interpretation might also be retained far from the classical limit, 
but then the Hamiltonian departs from the classical one, receiving quantum 
corrections. 

The extension to the QFT of a real scalar field <p in d + 1 (Euclidean) 
dimensions is quite straightforward. Let <p(x) = <p(r, x) where x = (r, x) e 
R( d+1) , rGK and x e R {d) . The free Euclidean action, So, defined in terms 
of the phase-space variables is 



Z*{P) = n B (uj) + 1 



1 



(30) 



1 




(31) 



with 



ft (7r,</?) = -[tt 2 + |V</f + mV 



(32) 
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We then have to implement the periodic boundary conditions both for 
<p(T, x) and its canonical momentum vr(r, x) 



<p(P,x) = <p(0,x) , vr(/5,x) = tt(0,x) , Vx eP, (33) 

which requires the introduction of two static Lagrange multiplier fields: £ a ( x ), 
a = 1, 2. Defining the two-component field <3> = ($ a ), a = 1, 2, such that 
$i = (p and $2 = 7r, an analogous procedure to the one followed for the 
harmonic oscillator yields: 

Zo(p) = u J vi Jv$ e -y dd+i x^ b ^+*f dd+i ^» , (34) 

where j a (x) = f a (x)(5(r - /3) - 5(r)) and: 

£ = ( 4 * ) ■ (35) 

where we have introduced h = \J — V 2 + m 2 , the first-quantized energy oper- 
ator for massive scalar particles. 

Performing the integral over yields the partition function in terms of 
the Lagrange multipliers: 

Zq(J3) = j ^ e -5/^'/ d<i ^ a (x)<x|M af) | y ) 5i) ( y ) ? (36) 

with M = 2(1(0) - £l(/3) - fl(-P) and 

)h- 1 fsgn(r) \ _ m 



Q(r) = 2 , , 2 °7 W e^ |T| . (37) 
Then, 

^ = ( V 1! ) ( ^ + irl ' (38) 

where 

n B = -J— . (39) 
Coming back to the expression for Zq(P), we see that: 

Z (f3) = Jvi^ V {- l -Jd d xJ d d y[U*) (xl^ 1 {n B + l)- 1 ^} fc(y) 
+ 6(x) (x|/ i (r iB + l)- 1 |y)6(y)]} • (40) 



By a simple field redefinition, we see that: 



Zo(l3) = det (n B + 1) (41) 
which can be evaluated in the basis of eigenstates of momentum to yield: 

MP) = II M^k) + !] (42) 



where = y/k 2 + m 2 . The free-energy density, Fq(/3), is of course: 

F °(« = js/S h,(1 " e "'* ) - (43) 

In the classical, high-temperature limit, the path integral for the partition 
function becomes: 

Z ((3) ~ Jvte-P*®, (44) 

where: 

#(0 = 5/ ^K(x) + |V6(x)| 2 + m 2 £ 2 2 (x)] . (45) 

This is, again, the usual classical expression for the partition function, with 
the Lagrange multipliers playing the role of phase space variables, and the 
integration measure being the corresponding Liouville measure. Besides, it is 
clear that the representation (l4"U|) always involves static fields, unlike in the 
Matusbara formalism. The price to pay for this 'dimensional reduction' is 
that the resulting 'action' (the exponent of the functional to be integrated) 
is spatially non local. It becomes local only in the high-temperature limit. 

When self-interactions are included, instead of the free action So, we must 
consider instead 

S = So + Sj , (46) 

where we assume that Si = f d d+1 xV(ip), V(cp) being an even polynomial in 
if. The new representation for the partition function is then: 

Z((3) = M j Vi J £>$ e-sm + ifd^xj^ _ ( 47 ) 

Recalling the definition of the generating functional for T = connected 
correlation functions, W, we may write: 



M I £>$ exp 



-S($) + % I d d+1 xj a § 



e~ w( - j) , (48) 
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so that the partition function Z(f3) becomes: 



Z{(3) = J Vi e~ w «) , (49) 

a formula that makes contact with T = QFT, building the thermal cor- 
rections on top of the T = ones. We have use the small j to denote the 
2-component current which is a function of the Lagrange multipliers. A capi- 
tal J shall be reserved to denote a completely arbitrary 2-component source, 
so that: 



J\f J T>§ exp -£($) + % J d d+1 xJ a $ a 



e~ w ^ . (50) 



Defining now H e ff(£), the 'effective Hamiltonian' for £, by means of the 
expression: 

Heff(0 = ^W(j), (51) 
we see that the partition function is simply given by: 

Z(P) = Jvi exp [-/?iZe//(0] > ( 52 ) 

one of the most important results in this article, since it yields the path 
integral for the quantum partition function as a classical-looking functional 
integral, involving an effective Hamiltonian which takes into account all the 
T = quantum effects. Indeed, the usual functional expansion for W(J) is: 

oo 

W(J) = Y± d^X^'-d^Xn 
n=2 J 

X WtLjXb ...,!„) J ai {xi) ■ ■ ■ J an {Xn), (53) 



where is i n times the n-point connected correlation function for the 

field $. Thus, knowing those correlation functions, at a given order in a loop 
expansion, one may obtain an analogous functional expansion for the effective 
Hamiltonian. Indeed, to do that, one should perform the integrations over 
the imaginary-time variables (taking advantage of the 5- functions): 

oo 

HefM) = E- /rf d X 1 ---^X n i/W an (x 1 , ...,X n ) 
n=2 ' 

x ^ 1 (xi)-.-^ 0n (x B ) , (54) 
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where 



a\...a r . 



x (5{t 1 -/3)-6(t 1 )) ... (S(r n - (3) - 8(r n )) . 



Of course, one will usually evaluate just a few of the terms in the ex- 
pansion for W, and as a consequence will obtain an approximation to the 
effective Hamiltonian. From the first non-trivial order, namely, keeping only 
the first term, the 2-point correlation function @ we see that: 



where 



and 



i J dx 1 dx 2 ^a{x.i)H^\x 1 , x 2 )£ 6 (x 2 ) 



Hj£(*i, x 2 ) 



1 r 



(2), 



- (J3, x a ; 0, x 2 ) - W% (0, Xl ; (3, x 2 ) 



(2), 



(55) 
(56) 



(57) 



Introducing Fourier transforms for the correlation functions W a b and the 
kernel, , ( !57l ) implies: 



^?(-,k) = i y_ ^[i-cos(M]w^(c,k). 



2^ 



(58) 



(2) 

On the other hand, the kernels W ah can be obtained from the T = 

(2) 

correlation function Wjy which is, essentially, the full quantum propagator 
for the </?-field: 

W{?(xi,x 2 ) = -Mxi)^(x 2 )> • (59) 

(2) 

Indeed, the explicit relations between the correlation functions W ah with 
a 7^ 1 or b ^ 1 and , in momentum space, are: 

W^^k) = 1 - c^wg^k) 
W{?(u;,k) = -wWjJ^k) 



(60) 



2 The meaning of this approximation shall be discussed at the end of this calculation. 
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In order to evaluate the frequency integrals, we need of course to make some 
assumptions about the structure of Wq • Assuming that there is only one 
stable particle, and that the renormalization conditions corresponding to the 
physical mass and wavefunction renormalization have been imposed. Then 
VVii has only two single poles, at the points ± i E ren (k), where E ren (k) is the 
dispersion relation, with a renormalized mass, and including all the quantum 
corrections. 

Considering each component of equation (1551 ) for specific values for a and 
6, we see that their respective to integrals are easily evaluated under the 
previous assumptions. The results are: 

(61) 
(62) 

while H$ = H% = 0. 

Since the functional integral has the same structure as the one for the 
free case, to evaluate Z(f3) in this approximation we just need to modify the 
expression for the energy in the free field result: Z{(3) ~ Z^i^P) where 

Z {2 \(3) = J] [n B (E ren (k)) + 1] (63) 

k 

which yields the free energy 

Fi2>i8) = \ J '^k Ml ~ e ~" E "" m) ■ m 

This expression looks like the free energy for a free field, with the renormal- 
ized energies instead of the free ones. This approximation amounts to consid- 
ering the particle states as non interacting, after the self-energy corrections 
have been taken into account (in the 2-point, T = correlation function). 
A possible way to understand this approximation is by the large- N limit, 
whereby correlation functions involving more than 2 points are suppressed. 

An interesting approximation scheme, related to the functional expansion, 
is given by the High-Temperature expansion. In terms of H e ff, we see that 
the only place where (3 appears is in the source j. Performing an expansion 
of j in powers of (3: 

00 on 

Ja(x) = £,(x) £(-l)»*LaW( r ), (65) 

n=0 



H[f(k) = E-i(k) n B (E ren (k)) 



Hg\k) = E ren (k) n B {E Ten (k)) + 1 
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and inserting this expansion into the expression for H e ff, we see that the 
kernels corresponding to its functional expansion become: 

flniH hrc m -l 

H (m) fxi ••• x ) = V - 

rai=l,— ,n m 

x [d^-.-d^W^Jx,,--- ,x m )] T ^ Q . (66) 
When /3 — > 0, the leading term is linear in (3 and only H^> survives: 

^) 2 ( Xl)X2 ) ~ /5 2 [9 Tl 9 T2 WS 2 (x 1 ,x 2 )] riiT ^ . (67) 
Keeping higher order N terms in the expansion: 

N 

Hi e?f(0 = E ~ { / ^ • • -^#£LjXl, • • • > X n ) 
n=2 ' J 

x ^ ai (xi)---^ n (x n ) , (68) 
the partition function naturally admits a perturbative expansion, with the 

(2) 

quadratic kernel H ab determining the propagator, and the higher order terms 
the (nonlocal) vertices. 

The final example we consider is a massive Dirac field in d + 1 space- 
time dimensions. The procedure will be essentially the same as for the real 
scalar field, once the relevant kinematical differences are taken into account. 
The action Sq for the free case is given by Sq = J d d+1 xip(ft + m)ip where 

# = ll = 7 M and {7 M , lv\ = 26^. 

We then impose antiperiodicity conditions for both fields: 

^(Ax) = -^(0,x) , ^(/3,x) = -^(0,x) (69) 

as constraints on the Grassmann fields. Those conditions lead to the intro- 
duction of the two S— functions: 

x exp [ - S£$,V0] • (70) 

Since the Dirac action is of the first-order, the introduction of two constraints, 
and two Lagrange multipliers, appears in an even more natural way than for 
the previous case. Those auxiliary fields, denoted by x( x ) and x(x) must be 
time-independent Grassmann spinors. The resulting expression for Zq{(3) is 
then 

Zl{0) = Af J V x VxV^Vije- s ^^ )+i ^ dd+lx ^ + ^\ (71) 
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where r] and 7] are (Grassmann) sources depending on \ and x through the 
relations: 

V (x) = X (x) [6(t -f3) + 5(r)] , f}(x) = *(x) [6(r - 0) + 5(r)] . (72) 
Integrating out ip,4>, we arrive to: 



VxDx exp 



P H *ff(x,x) 



where 



with: 



H eff (x,x) = j d " x j Ax(x)# (2) (x,y)x(y) 



(73) 
(74) 



^ 2 )(x,y) = (x,0|(^ + m)- 1 |y,0) + (x,/5|(^ + m)- 1 |y,/5) 
+ (x,0|(^ + m)- 1 |y,/5) + (x,/5|(^ + m)- 1 |y,0) 

= - ^/(O.x-y) +5/ (/3,x-y)+5 / (-/3,x-y)](75) 

In the last line, S 1 /, denotes the Dirac propagator. A quite straightforward 
calculation shows that 



1 



(76) 



where h F = ( 1 + e^) is the Fermi-Dirac distribution function, written 



in terms of with h, the energy operator (defined identically to its real scalar 
field counterpart), u is a unitary operator, defined as 



h 



a 



D 



h D ee 7 • V + m . 



Then we verify that: 

Z*(p) = detu det _1 [(l -h F )l] , 
(I ee identity matrix in the representation of Dirac's algebra) 



(77) 



(78) 



z£(P) = \ U\ 1 + e 



-(3E(P) 



(79) 



with E(p) = \Jy> 2 + m? and ee dimension of the representation (we have 
used the fact that det u — 1). 
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Again, the procedure has produced the right result for the partition func- 
tion, with a normal-ordered Hamiltonian. On the other hand, for a Dirac 
field in a static external background corresponding to a minimally coupled 
Abelian gauge field the A = gauge, we have 



d d+i 



x 



i^(x)( $ + ie'j ■ A(x) + m)ijj(x) 



(80) 



The assumed r— independence allows us to carry on the derivation described 
for the free case, with minor changes, arriving to the expression: 



Z f (f3) = detu(A)det- 1 (n F (A)l) 

(l + e-W)l 



e^ (A) det 



where 
and: 



h(A) = V-D 2 + m 2 , D = V — ieA 
det (7 • D + mj 



JK(A) 



(81) 



(82) 



(83) 



Notice that the factor det 



det \J — D 2 + m 2 

can be formally diagonalized in 



(l+e-« A ')l 

terms of the energies E\(A) in the presence of the external field. Thus we 
arrive to the expression: 



e iK(A) x 




-I3E X (A) 



(84) 



The factor e tK ^ A \ on the other hand, is topological in origin, as it depends 
on the phase, K(A), of the determinant of ho- On the other hand, he 
may be regarded as a kinetic operator in one fewer dimension. For Dirac 
fermions, we know that the phase of det ho can be non-trivial only when d 
is odd, i.e., when d+ 1 is even. However, the 7-matrices appearing in det^ 
form a reducible representation of the Dirac algebra in d dimensions, with 
the matrix 7 T relating every eigenvalue to its complex conjugate. Thus, as a 
result, the phase K(A) vanishes. Of course, a non-vanishing result may be 
obtained for other fermionic systems, like Weyl fermions for d + 1 = even. 

Summarizing, we have shown that, by introducing the 'thermal' boundary 
conditions as constraints in the Euclidean path integral for the vacuum func- 
tional, we may obtain a novel representation for the partition function. This 
representation may be thought of as an integral over the phase space vari- 
ables, weighted by a Boltzmann factor corresponding to an effective quantum 
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Hamiltonian, H e jf, which reduces to the classical one in the corresponding 
(high-temperature) limit. We discussed the main properties of this represen- 
tation for the cases of the real scalar and Dirac fields, two typical examples 
that have been chosen for the sake of simplicity. It is not difficult to general- 
ize the representation to the case of systems containing fermions interacting 
with bosons, by introducing the proper T = generating functionals. 
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